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We show that negative of the number of floppy modes behaves as a free energy for both con- 
nectivity and rigidity percolation, and we illustrate this result using Bethe lattices. The rigidity 
transition on Bethe lattices is found to be first order at a bond concentration close to that predicted 
by Maxwell constraint counting. We calculate the probability of a bond being on the infinite cluster 
and also on the overconstrained part of the infinite cluster, and show how a specific heat can be 
defined as the second derivative of the free energy. We demonstrate that the Bethe lattice solution 
is equivalent to that of the random bond model, where points are joined randomly (with equal prob- 
ability at all length scales) to have a given coordination, and then subsequently bonds are randomly 
removed. 



I. INTRODUCTION 

Connectivity percolation on Bethe lattices or infinite 
statistically homogeneous Cayley trees was thoroughly 
analyzed by Fisher and Essam, [Q and provides a useful 
model for percolation. Rigidity on Cayley trees was first 
studied by Moukarzel, Duxbury and Leath; ||] henceforth 
referred to as MDL. This model has g degrees of freedom 
per site, with g = 1 corresponding to connectivity per- 
colation being a special case, and so may be regarded as 
an extension of the work of Fisher and Essam. In this 
paper, we develop a free energy for the rigidity percola- 
tion problem, and as an example of its use, show how to 
locate the bulk rigidity transition on the Bethe lattice. 
We use the term Bethe lattice, rather than Cayley tree, 
to emphasize the bulk behavior of lattices containing no 
loops away from the boundary. ||,p[ 




FIG. 1. Showing a tree with a single bar b — 1 along 
the 2 = 3 bonds from each site. The bonds are present with 
probability p. Each site, shown as a shaded circle, has g 
degrees of freedom. Only the first few levels from the rigid 
busbar are shown. 

We have previously suggested that a free energy can be 
defined as the negative of the number of floppy modes. 

In this paper we prove this assertion for the case of 
random bond dilution in a general lattice with g degrees 
of freedom per site, and with z nearest neighbors, con- 
nected by 6 = 1 bars. For g = 1, this gives connectivity 
percolation. For g = d, this corresponds to the central 
force model [^,^, where g = 2 \n two dimensions and 
g ~ Q in three dimensions. For bodies joined by bars, 
g = d{d -I- l)/2 in a d dimensional space. ||] 

We prove that this free energy has the correct convex- 
ity property and hence a specific heat can be defined from 
the second derivative. This definition of the free energy 
reduces to the known result for connectivity percolation 



when f; = 1, where it becomes the negative of the total 
number of isolated clusters. This latter result can ei- 
ther be found directly, j^] or as the s ^ 1 limit of the s 
state Potts model. Q For rigidity, this free energy allows 
us to locate the bulk transition on the Cayley tree, and we 
refer to this as the Bethe lattice solution. [|| We confirm 
that for b = 1 and g > 2 the rigidity transition is always 
first order on the Bethe lattice, whereas of course it is 
second order for connectivity percolation where g ~ 1. 
On the triangular lattice, rigidity percolation is second 
order. whereas it is first order on the Bethe lattice. 

The nature of the rigidity transition is a subtle question, 
depending on the network, and this has not always been 
taken into account, p^ ] For connectivity percolation, the 
transition is always second order. ||l|,[l3| 

MDL developed a general solution, using a transfer ma- 
trix technique, for a network of rigid bodies, each with 
g degrees of freedom generically connected by b bars, as 
sketched in Fig. 1. Initially each rigid body has z neigh- 
bors, and the rigidity of the network is studied as the 
bonds are randomly removed. For b < g, this type of 
network is always floppy, unless the tree is attached to a 
rigid surface (busbar) , in which case rigidity may or may 
not propagate away from the busbar, depending on the 
degree of dilution. 

We amplify the work of MDL by using the free energy 
and the associated Maxwell construction, Q to remove 
boundary constraints in the analysis of bond diluted Cay- 
ley trees, and hence obtain the bulk or Bethe lattice so- 
lution. This constraint equation leads to a higher rigidity 
threshold than that found in MDL. In analogy with ther- 
mally driven thermodynamic transitions, we interpret the 
threshold found in MDL as a spinodal point and use ps 
when referring to it in this paper. We use pc to refer to 
the bulk threshold. 

The layout of this paper is as follows. In the next 
section (Sec. II) we develop results applicable to general 
bond-diluted lattices, namely: the proof that the nega- 
tive of the number of floppy modes acts as a free energy 
for studying bond-diluted connectivity and rigidity per- 
colation, and we derive the associated Maxwell construc- 
tion or consistency condition. We review the main 
equations from MDL for the Cayley tree, and evaluate the 
probability Poo of being in the infinite rigid cluster and 
also the probability Pov of being in the overconstrained 
part of the infinite cluster; both defined in the bulk far 
from the busbar. 

In Sec. Ill, we discuss the random bond model, 
where z bonds join each site with other sites, regardless 
of distance, and then these bonds are randomly diluted. 
This model is solved using the integer algorithm, the peb- 
ble game, [|9|, [To| , p^ for the connectivity case g = 1, z — 3 
and the rigidity case g = 2, z = 6. We demonstrate that 
in the thermodynamic limit, the random bond model is 
equivalent to a Bethe lattice. 

In Sec. IV, we show that for g — 1, the familiar con- 
nectivity case is recovered. This analysis is examined for 
a general z, specializing to 2; = 3 to illustrate the detailed 



behavior for connectivity percolation. 

In Sec. V, we give detailed results for a Bethe lattice 
with g = 2 and z = 6 as an example of a first-order 
rigidity transition. We apply the consistency condition 
to locate the first order transition in the bulk, and show 
that the results are equivalent to the random bond model. 
We also give results for the probability Poo of being on 
the infinite cluster and the probability Po„ of being on 
its overconstrained part. We show that the results for 
the number of floppy are equivalent to the g — 2 random 
bond model. We present a table summarizing important 
results for a range of values of g and z. 

In Sec. VI, we show how the limit z — > cx3 can be taken. 
We show when a further limit — > cx) is also taken, we 
recover the Maxwell constraint counting result, which we 
therefore identified as the mean field theory for rigidity. 

Throughout this paper we focus our attention on bond 
dilution where the number of bars 6 = 1. Results for 
b > 2 and/or for site dilution can probably be obtained 
in a similar manner. 



II. FORMALISM 
A. General 

In this subsection we develop a free energy for both 
connectivity and rigidity problems on general lattices, 
and show that this leads to a consistency constraint on 
the number of floppy modes that is useful in locating first 
order transitions. These results are applied to Cayley 
trees in the next sections. 

A constraint counting method (also called Maxwell 
counting) has been very useful as a conceptual tool in 
understanding rigidity, and in particular the rigidity of 
glasses. P,pj|-p^ The Maxwell counting approach centers 
on the number of floppy modes, where we define f{p) as 
the average number of floppy modes per degree of free- 
dom. We consider generic networks where each site has g 
degrees of freedom with central forces that connect to z 
neighbors. Bonds are randomly present with probability 
p. An exact relation on any bond-diluted lattice is 
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where F is the total number of floppy modes. There 
are N sites in the network, and r{jp) is the probabil- 
ity that a bond is redundant, i.e. its removal does not 
change the number of floppy modes in the system. Eq. (|^) 
demonstrates that we can count either floppy modes or 
redundant bonds, and although early work has empha- 
sized counting floppy modes, 18 1 they are equivalent. 
In Maxwell counting, r{p) is assumed to be zero for 
p < Pm and the rigidity threshold is taken to occur at 
f{Pm) — 0, which leads to the Maxwell estimate of the 
rigidity threshold. 



Pm. = '2g/z. 
Within this approximation it is evident that 
f{p) = for p>Pm 

f{p) = 1 - p/p„i for p < Pra, 



(2) 



(3) 



so that f{p) changes slope at p,„. We emphasize that 
Maxwell counting is not correct, but does provide a useful 
initial approximation. In reality f =/= for p > p„i and 
is non-zero right up to p = 1 as there is always a small 
probability of having a floppy inclusion if p 7^ 1 . This is 
a Lifshitz type of argument. |20 

We note that /(O) = 1 and /(I) = 0, and hence from 
Eq. (H) we have r(0) — and 



r(l) = l 



2.9 



(4) 



Equivalently, the following relation must be satisfied 



/«dp = -l. 



(5) 



These consistency conditions remove boundary con- 
straints from the calculation and hence enable prediction 
of bulk critical behavior. This is a key new element of the 
Cayley tree theory, which was lacking in MDL, and which 
we now refer to as the bulk or Bethe lattice solution. In 
particular it allows us to locate the bulk critical point 
Pc when the transition is first order. In MDL, the point 
at which the finite real solution to the mean field equa- 
tions ceased to exist was identified as the critical point. 
Here we label that point as Ps and reinterpret it as the 
spinodal point in analog with thermally driven first-order 
transitions. 

We now show that the number of floppy modes acts as 
a free energy in both connectivity and rigidity problems. 
To demonstrate this, we recall the relation developed by 
Jacobs and Thorpe [|ol 
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where / is the number of floppy modes per degree of 
freedom, Nq is the number of overconstrained bonds on 
the lattice, Nj is the number of isostatic bonds on the 
lattice, and 
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is the total number of bonds present on the lattice, and 
the first derivative /*^^^ — df /dp. The relation is 
obtained by removing a single bond, chosen at random, 
and ascertaining the probability that this bond is in 
an overconstrained region (in which case the number of 
floppy modes is unchanged) or else the number of floppy 
modes is reduced by 1. That is (for one removed bond) 



AF = (1 — Nq/Nb), and if we now remove ANb bonds, 
we have 



B. Cay ley Trees 
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and the result follows. By comparing Eqs. (|i|) and 
, we see that the number of overconstrained bonds is 
given by the rate of change of the number of redundant 
bonds. 



drjp) ^ No 
dp Nb 



(9) 



We can use a similar argument to derive the important 
result 
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where the second derivative /^^^ = f /d^p. We demon- 
strate this as follows. If we remove one bond, chosen 
at random, then the change in the number of overcon- 
strained bonds is given by 



AA^o = -^(l + A) 



(11) 



where A > because when an overconstrained bond is 
removed, the number of overconstrained bonds is reduced 
by at least one, and sometimes by more. If we remove 
ANb bonds, we therefore have 
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ONb Nb 

Using this, we find 
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(12) 
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and hence using (9) and (10), 
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establishing that —f{p) is a convex function of the frac- 
tion of bonds present p. Note that we must ensemble 
average the results (^2|) and (14) so that A is to be inter- 
preted as an ensemble averaged quantity. 

Therefore we can use —f{p) as a free energy, and if 
there is any ambiguity the system will always be in the 
lowest free energy (maximum floppy modes) state. The 
quantity —f(j>) is convex as required of a free energy, 
and we will refer to f{p) interchangeably as the fraction 
of floppy modes or as the free energy. Note that the 
proof above is given for arbitrary g and z, but restricted 
to 6 = 1 and bond dilution. Because f^'^^{p) is the second 
derivative of a free energy and is positive definite, we can 
regard it as a specific heat and it is calculated explicitly 
in subsequent sections, pl^] 



Following MDL, we consider trees which have co- 
ordination number z with g degrees of freedom per site, 
which form a Cayley tree network attached to a rigid 
boundary which we call a busbar. This is shown in Fig. 1. 
The busbar is not necessary for connectivity when .g = 1 , 
but is of vital importance for rigidity when g > 2. We 
define Tq to be the probability that a bond on a branch 
n levels away from the busbar is part of the infinite rigid 
cluster. In general, if the sites of the tree have g degrees 
of freedom, rigidity is transmitted to the next level of the 
tree provided at least g of the bonds to the lower level 
are occupied and provided that the sites at the ends of 
these bonds are rigid. This gives the recurrence relation 
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where Tq is the probability a bond n levels from the 
busbar is rigid. If we take the thermodynamic limit 
(very large n) , Eq. ( p^ ) iterates to a steady-state solution, 
which we call Tq and is given by 
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From this equation, we can find the probability of having 
a single degree of freedom with respect to the (distant) 
boundary 

Ti^(^~_^^{pToy-\l-pToy-', gyi. (17) 

and more generally for / degrees of freedom with respect 
to the boundary 
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Summing over all possibilities, we have the useful sum 
rule 
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Equation ( p^ ) is the self-consistent equation for the rigid- 
ity order parameter on bond-diluted Cayley trees. MDL 
also considered a more general class of problem in which 
another degree of freedom & > 2 [the number of con- 
straints (or bars) between each pair of sites] is allowed. 
They also considered the case of site dilution, which is 
trivially related to bond dilution on a Cayley tree for 
h — 1 only. 

For trees with & = 1, the probability that a bond is 
overconstrained is pTg , i.e. a bond is overconstrained if 



it is present and both of the sites at its ends are already 
rigidly connected to the busbar. Such bonds are the 
overconstrained bonds of the infinite cluster, as no other 
bonds are overconstrained far from the busbar, which is 
the bulk solution that we are seeking. Thus from we 
have, 
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Differentiating a second time gives 

9 dp 

Since To > 0, we must have solutions where 



dp 
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to be acceptable. 

In a similar way the probability that a bond is con- 
nected to the infinite rigid cluster via the busbar at one 
end, but has one degree of freedom at the other is 2pToTi. 
Such bonds are the isostatic bonds associated with the 
infinite cluster. p|-pl] Thus the total probability that a 
bond, that is present, belongs to the infinite cluster Poo 
is 



Poo^To' + 2ToTi. 



(23) 



The probability that a bond that is present and belongs 
to the overconstrained part of the backbone is given by 



Pr, 



To' 



(24) 



Since the number of redundant bonds r(p) is zero for 
p < Pc and after finding r^^^ (p) = Tq from (20) using (1), 
we must have, 
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In particular, combining Eq. (25), with the consistency 
condition (4), we have 
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(26) 



This is the key new component of the theory developed 
here. As in MDL, we solve the tree equations for the 
order parameter Tq but now we use Eq. (26) to identify 
the rigidity threshold Pc- Once we have Pc, the number 
of floppy modes at any p is given by [from Eq. (1) using 
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The condition, Eq. (26) may be regarded as a Maxwell 
condition, in analogy with that commonly used to locate 
the first order transition in thermodynamic systems, p^ ] 



III. RANDOM BOND MODEL 

The random bond model consists of N sites thrown 
down randomly, with each site joined at random to z 
other sites (without regard to distance) to form a z- 
coordinated network. Small rings appear with proba- 
bility 0(1/A^). It is convenient to show this network in 
a plane as in Fig. 2, where z = 3. The dimension of this 
model is only defined through the number of degrees of 
freedom g associated with each site. In the Bethe lattice 
and in the random bond model, there are no loops in the 
thermodynamic limit N oo and so in a diagrammatic 
expansion, all terms agree, and the models are equiva- 
lent. For finite Cayley trees, there are loops involving 
the busbar as can be seen in Fig. 1. Likewise, there are 
loops in the random bond model as can be seen in Fig. 2. 
So the equivalence is only in the bulk thermodynamic 
limit. We find this is a useful alternative viewpoint on 
the Bethe lattice limit. 

We have numerically examined connectivity percola- 
tion having g = 1, z — and the g = 2, z = 6 rigidity 
percolation case using the random bond model. These 
simulation results are compared with the corresponding 
exact Bethe lattice calculations in the next two sections. 
The pebble game ^,|l^,|l^ was used to find the number 
of floppy modes and the derivatives /'-^-'(p) and /^^Hp) 
for the g — 2, z — 6 case. For connectivity percolation 
with g = 1, a similar pebble game algorithm can be con- 
structed where only one pebble is assigned to a site. |2|] 
In fact, as long as a site represents a rigid body having g 
degrees of freedom and not a point, the pebble game can 
be straight forwardly generalized by assigning g pebbles 
to each site. 23 

For the connectivity percolation case, however, we ac- 
tually used the g = 2 pebble game, (two pebbles per site), 
by invoking an interesting mapping that makes g = 2 
rigidity percolation equivalent to connectivity percola- 
tion. This mapping is valid for any generic network and 
is not limited to the random bond model. The mapping 
consists of adding a single ghost site to the N site net- 
work. The ghost site has = 2 as do all other sites in 
the network. An additional N bonds are placed between 
each of the N sites in the network and the ghost site. 
With this additional ghost site and its associated bonds, 
the network consists of edge sharing triangles. That is, 
any three sites in the network that are connected, are 
mutually rigid because they form edge sharing triangles 
with the ghost site common to all triangles. 

The ghost site and its associated bonds, allow connec- 
tivity to be a sufficient condition for a set of sites to be 
mutually rigid. Since rigid clusters are those that are 
connected, a one to one mapping between clusters in a 
g = 2 network to a g = 1 network is established. How- 
ever, the number of floppy modes of the g = 2 network 
(with the ghost site) will be two more than the num- 
ber of rigid clusters. These floppy modes can be viewed 
as two translational motions for the ghost site, and one 



rotational motion of each rigid cluster about the ghost 
site. By accounting for these two extra modes, an exact 
mapping is established. 




FIG. 2. A sketch of the random bond model with a single 
bar 6=1 along the z = 3 bonds from each site. The bonds 
are present with probability p. Each site, shown as a shaded 
circle, has g degrees of freedom. Only A'^ = 12 sites are shown 
in this sketch. In simulations using the pebble game, we use 
iV = 262, 144. 



on the infinite cluster is given from Eqs. (|2 
and using Ti = 1 — Tq 
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When p = 1, we have To = 1, so that the ratio ( |3l| ) 
becomes unity, meaning that the entire tree is overcon- 
strained (no dangling ends). Near p = pc, and using 
( p9| ) , the fraction of overconstrained bonds in the infinite 
cluster goes to zero linearly as 
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so that as the connectivity percolation transition is ap- 
proached, the infinite cluster becomes completely iso- 
static. For z = 3 we can solve the quadratic equation 
(H) to get 



Tq ~ for p < p, 
{2p-l 



Tn = 



for p> Pc 
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where Pc = 1/2, and is shown in Fig. 3. 



IV. CONNECTIVITY PERCOLATION 



0.6 



The analysis in Sec. II holds for any value of g. For 
g = 1 we recover and extend the familiar results of con- 
nectivity percolation ||l|,|l3| for random bond dilution on a 
tree with coordination z. In this section we use a general 
z as much as possible, but then focus on z = 3 to simplify 
the algebra when necessary. For g = 1, Eqn. (|l^) reduces 
to 
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Expanding near pc it is easy to show that the order pa- 
rameter To vanishes as 
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The number of redundant bonds r(p) can be found from 
(^) by using (^sj) to first calculate dTo/dp and then 
doing the integration to give 



ip) = To 
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where To is obtained by solving the polynomial (|28| ). No- 
tice that the fraction of bonds in overconstrained regions 
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FIG. 3. The result for bond diluted connectivity percola- 
tion on a Bethe lattice for To is shown for co-ordination num- 
ber z = 3 and with g = 1 degrees of freedom per site. The 
probability that a bond is present is p, and percolation occurs 
at Pc = 0.5. 

We can use (^^ to check that we get the same value of 
Pc, by noting that we correctly reproduce r(l) = 1/3 [see 
Eq. (^]. Therefore the transition is indeed second order. 
Having found pc, we find f(p) and r{p) from Eqs. ( [25| ) 
and (I27I). 
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and differentiating again 
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FIG. 4. Results for bond diluted connectivity percolation 
on a Bethe lattice for the number of floppy modes f{p) and 
the first two derivatives. The number of redundant bonds 
r(p) are also shown. The second derivative f^^\p) acts as a 
specific heat for this problem. Results are for co-ordination 
number z = 3 and with g — 1 degrees of freedom per site. The 
probability that a bond is present is p, and percolation occurs 
at Pc = 0.5. The open circles are from computer simulations 
of the random bond model using the pebble game as described 
in Sec. Ill, using A'' = 262, 144 sites and averaging over 2,000 
realizations. The Maxwell estimate [Eqs. (2) and (3)]for the 
number of floppy modes is a straight line follows the f{p) at 
small p. 



These results ( p^ ) - (^^ are shown in Fig. 4. We note 
that this second derivative, which is like a specific heat, 
is strongly peaked around p — I — l/-\/6 = 0.59, although 
it does go zero at pc — 0.50. Note also that /, /^^^ and 
/*^^) are all continuous at pc, with /^'^^ being the first 
derivative to show a discontinuity. 

The number of overconstrained bonds in the infinite 



cluster is Pov and is given from (24) by 
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given by (g3|) 

Poo = 1 
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for p>Pc 



(39) 



which has been derived by Essam and Fisher, |]l|[see their 
Eq. (35b)] although notice they have an extra factor of p 
as they normalized to all bonds, not just those present. 
These results are shown in Fig. 5. The reader might be 
concerned that some overconstrained bonds are not asso- 
ciated with the infinite cluster. This is not so on a Cayley 
tree in the asymptotic limit where only the infinite cluster 
reaches back to the busbar. All other clusters are isolated 
and contain no overconstrained bonds, except finite rigid 
clusters attached to the busbar which are considered to 
be surface effects and hence are ignored in the bulk or 
Bethe lattice solution. 
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FIG. 5. Results for bond diluted connectivity percolation 
on a Bethe lattice for the probability of being on the infinite 
cluster Poo and the probability of being on the overconstrained 
part of the infinite cluster Pov Results are for co-ordination 
number z — 3 and with g = 1 degrees of freedom per site. 
The probability that a bond is present is p, and percolation 
occurs at pc = 0.5. 



In the connectivity case f^^\p) and the order param- 
eters To, Pov and Poo are continuous at the transition, 
and so connectivity percolation is of a conventional sec- 
ond order type, as indeed it is on real lattices like the 
triangular network. p3[ | 

V. RIGIDITY PERCOLATION 

In the rigidity case we must numerically solve Eq. ( p^ ) 
for To and hence find the other quantities of interest. We 
first solve the self-consistent Eq. (|l^), by simple itera- 
tion, to find To, which is shown in Fig. 6 for g = 2 and 
z — 6. The analogue of the Maxwell construction (^) 
is then used to find Pc- This is shown as the jump in 
Fig. 6 at Pc = 0.656. The other possible solutions are 
ruled out by application of (^) and hence we have a first 
order jump from a rigid to a floppy state. The point at 
Ps — 0.603 can be interpreted as a spinodal point. ||2j] 
The dashed line is unstable as can be seen from Eq. (22) . 
Having found Tq and Pc, we can find f{p), /'^^ (p), /^^^ (p) 
and r{p) as shown in Fig. 7. It can be seen that there 
are no redundant bonds for p < Pc and hence f^^^ and 
/(^) are fiat in this region, as also happens for connec- 
tivity percolation as shown in Sec. IV. In Fig. 7, we also 
show the results from the random bond model, discussed 
in Sec. Ill, and find that the results are equivalent as 
expected. Note that with the random bond model, the 
first order transition occurs naturally and no Maxwell 
construction is needed, as the pebble game always finds 
the bulk equilibrium solution. 
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FIG. 6. Results for bond diluted rigidity percolation on a 
Bethe lattice for To are shown for co-ordination number z — 6 
with g = 2 degrees of freedom per site. The probability that 
a bond is present is p and percolation occurs at pc = 0.656 
shown by the vertical line. The thin line extends out to the 
spinodal point at pa = 0.603 and the dashed line shows the 
unstable solution. 



In Fig. 8, we show results for Poo and Pov from 
Eqs. (23)and (24) respectively. It can be seen that most 
of the infinite cluster is overconstrained, even at the first 
order transition. 

In the connectivity case as shown in Figs. 3 and 5, 
the order parameter [ To, Poo or Pov] is continuous at 
Pc, whereas in the rigidity case as shown in Figs. 6 and 
8 the order parameter has a large first-order jump at 
Pc- The quantity f'^^^ip), which acts like an energy, also 
has a large first order jump. One could argue that in 
the rigidity metastable rigid state exists for any 

Pc > p > Ps- Note that for connectivity percolation, the 
boundary conditions are irrelevant, and no difference is 
found with or without a rigid busbar. This is not so for 
rigidity percolation, where there are three real solutions 
p5| for all Ps < p < 1, with an unstable solution at 
finite To existing in this regime as shown by the dashed 
line in Fig. 6. In order for the stable finite To solution 
to be found, the boundary rigidity must lie above the 
unstable fixed point. Any boundary rigidity below this 
value iterates to the To = solution. This curious effect 
of the boundary condition in rigidity percolation, even 
on the bulk Bethe lattice solution, needs further study. 
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FIG. 7. Results for bond diluted rigidity percolation on a 
Bethe lattice for the number of floppy modes f{p) and the first 
two derivatives. The number of redundant bonds r{p) are also 
shown. The second derivative Z'^-* (p) acts as a specific heat for 
this problem. Results are for co-ordination number z = 6 and 
with g — 2 degrees of freedom per site. The probability that 
a bond is present is p and percolation occurs at pc = 0.656. 
The open circles are from computer simulations of the random 
bond model using the pebble game described in Sec. Ill, using 
A'' = 262, 144 sites and averaging over 2,000 realizations. The 
Maxwell estimate [Eqs. (2) and (3)]for the number of floppy 
modes is a straight line follows the f{p) at small p. 
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FIG. 8. Results for bond diluted rigidity percolation on 
a Bethe lattice for the probability of being on the infinite 
cluster Poo and the probability of being on the overconstrained 
part of the infinite cluster Pov Results are for co-ordination 
number z = 6 with <? = 2 degrees of freedom per site. The 
probability that a bond is present is p, and percolation occurs 
at pc — 0.656. 

It is useful to compare the results for pc with that 
predicted by Maxwell counting. Maxwell counting would 
predict pm = 2g/z = 2/3 in both cases studied here. 
In the connectivity case this is markedly different than 
the exact result pc = 1/2 because there are many floppy 
modes at pc in that case as can be seen from Fig. 4. 
In contrast, in the rigidity case, the Maxwell estimate 
is close to the Bethe lattice result Pc — 0.655 and there 
are rather few floppy modes at Pc as shown in Fig. 7. 
It is remarkable that Maxwell counting does so well in 
locating the position of both first order and second order 
P,p^ rigidity transitions. 

Results of calculations for a range of values of g and z 
for bond diluted Bethe lattices are summarized in table 
1. In all rigidity cases (i.e. g > 1) there is a large jump in 
the order parameter Tq, Pqo or Pqv at the critical point 
Pc, although the jump does decrease with increasing z for 
a given g. On real lattices, there will also be a Lifshitz 
tail [ pO[ in the number of redundant bonds r{p) [see 
Figs. (4) and (7)] extending into the regime p < Pc all 
the way to p = 0. Notice that the Bethe lattice result 
does correctly reproduce the Lifshitz tail in f{p) in the 
region < P ^ 1 • 

It is seen from Fig. 7 that the random-bond model is 
numerically equivalent to the Bethe lattice results, as ex- 
pected from our previous discussion. The pebble game 
is an exact algorithm for a particular network and gives 
the exact number of floppy modes (and the associated 
derivatives). There is therefore no hysteresis in the re- 
sults of the pebble game on the random bond model, and 



in these large samples, the system flips from the rigid 
state to the floppy state around pc — 0.656 as a single 
bond is removed in a given sample. This is a very un- 
usual situation as most numerical techniques, like Monte 
Carlo encounter hysteresis near a first order transition. 

VI. LARGE Z LIMIT 

It is instructive to take the limit as the number of 
neighbors z increases to infinity. In order to make this 
limit sensible, it is necessary to change from the p variable 
to r = zp, which is the mean coordination of a site. For 
connectivity percolation, with g = 1, Eq. (16) becomes 

To = 1 - e"''^" (40) 

which leads to a rc = zpc = 1 which differs from the 
Maxwell estimate of r„i — zpm = 2 by a factor 2. Of 
more intertest to us is the case g — 2 corresponding to 
rigidity, where Eq. (16) becomes 

ro = l-(l + rro)e-'-^°. (41) 

Eq. ( pl[ ) can be solved for To and the spinodal point 
found explicitly by differentiation. In this large z limit, 
the analog of the Maxwell condition given in Eq. (26) 
becomes 

/>oo 

/ {l^T^)dr = 2g~rc. (42) 

Eqs. ( |4l| ) and (^) can be solved simultaneously to locate 
the flrst order transition at and the result is listed 
in Table 1. We also show the results for the number 
of floppy modes / at the transition and at the spinodal 
point, using Eq. (27). The values of Pa^ and Pov are 
obtained using Eqs. (23) and (24) respectively. It can be 
seen from Table 1 that the values of all these quantities 
evolve smoothly as z — > c». 

For larger values of g, we use Eq. (16) to generalise 
Eq. (^ ) to arbitary g 

^0 = 1- (E^)^-^""- (43) 
\fe=o ■ / 

Notice that Eq. (||) reduces to Eqs. (|^ and (||) re- 
spectively for .9=1 and g — 2. Solving Eqs. ( [431 ) and 
( p^ ) simultaineously, we flnd the results for 5 = 3 and 6 
listed in Table 1. 

By examining the equations in this section, we can find 
that as g increases, the values of the mean coordination 
at the spinodal point Vs and at the first order transition 
Tc so that in the limit g = oo, we have 

2r, = ra = rrn = 2g (44) 

and we recover the Maxwell result given in Eqs. (2) and 
(3). Here is defined by r„i — zpm- It is interesting 



to note that the spinodal point is always distinct from 
the transition, and differs by a factor 2 in this hmit. We 
obtained the hmiting value of rs by carefully tracking Eq. 
(43) for large g, and find that the limit is approached very 
slowly. In this limit the metastable region between the 
spinodal and the first order transition exists between Ts 
and Tc, and we have the Maxwell result which is also first 
order of course. 

Note that in this limit g — > cxd, we have To = for 
r < Tc and Tq = I for r > Tc, so that the system is 
either perfectly floppy for r < Tc when there are no re- 
dundant bonds (this is true for arbitary z and g) and now 
also perfectly rigid for r > 7'c in this limit with z ^ oo 
and g ^ oo when there are are no floppy regions. This 
limit in which the system behaves homogeneously and 
the fluctuations are suppressed, may be regarded as the 
mean field limit for rigidity. 

VII. CONCLUSIONS 

In this paper we have shown that the free energy is 
negative of the number of floppy modes. This general- 
izes a previous result which gives the free energy as 
the total number of clusters in connectivity percolation, 
where there is a single floppy mode associated with each 
isolated cluster. In connectivity percolation, the free en- 
ergy can also be found as the limit as s — > 1 of the s 
state Potts model. |8| No equivalent approach has been 
possible for rigidity. 

We have used the free energy, and the associated 
Maxwell type construction, or self consistency condition, 
to show that the rigidity transition is always flrst order on 
Bethe lattices. It is necessary to use this construction in 
order to remove boundary constraints. In previous work, 
P,p5[ a boundary dependent transition was found at the 
spinodal point ps. As is expected from experience with 
thermodynamic transitions, metastability and boundary 
effects can be important p4| in the regime ps < p < Pc, 
and this needs further investigation. 

While rigidity concepts have been extensively applied 
to experiments in chalcogenide glasses first or- 

der jumps in experimental quantities have never been 
observed (although see recent Raman results in chalco- 
genide glasses p7|). 
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Table 1 This table gives a summary of results for bond-diluted Bethe lattices with g degrees of freedom per site 
which arc joined by 6 = 1 bars to z neighbors. The Maxwell estimate for the transition r„j = zp„i is given as well as 
the spinodal point = zps and the critical point rc = zpc, where we mulitply these quantities by the coordination 
number z to make them coordination numbers, that are semi- invariant. The number of floppy modes at the spinodal 

point is f{ps) and at the critical point is f{pc)- The jumps in the order parameters arc given for Tq, and the probability 
Poo that a bond is present and is part of the infinite cluster, and the probability Pov that a bond is present and is 
overconstrained. 



